Correlation induced spin freezing transition in FeSe: a dynamical mean field study 
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The effect of local Coulomb interactions on the electronic properties of FeSe is explored within 
dynamical mean field theory combined with finite-temperature exact diagonalization. The low- 
energy scattering rate is shown to exhibit non-Fermi-liquid behavior caused by the formation of 
local moments. Fermi-liquid properties are restored at large electron doping. In contrast, FeAsLaO 
is shown to be located on the Fermi-liquid side of this spin freezing transition. 

PACS. 71. 20. Be Transition metals and alloys - 71.27-fa Strongly correlated electron systems 



The recent discovery of high-temperature supercon- 
ductivity in iron-based pnictides ,lj] and chalcogenides ^] 
has led to an intense discussion concerning the role of 
Coulomb correlations in these materials. Although com- 
pounds such as FeAsLaO (1111), BaFe2As2 (122), LiFeAs 
(111), and FeSe (11) all have rather similar one-electron 
properties, a variety of experiments suggest significant 
differences. For instance, photoemission measurements 
show that FeAsLaO is moderately correlated, with about 
50 % Sd band narrowing and effective mass enhancement 
of about 2 to 3 3]. In contrast, several photoemission 
data on FeScajTei-^: samples reveal larger effective mass 
enhancement and stronger band narrowing 0^3' op- 
tical measurements exhibit incoherent spectral features 
indicative of a pseudogap , and several transport mea- 
surements show large deviations from Fermi-liquid be- 
havior 13" 13 1 . Moreover, recent theoretical work fT3| 



suggests that Coulomb interactions in (11) compounds 
ought to be less well screened than in (1111) systems. 

The aim of this paper is to elucidate the origin of the 
experimentally observed bad-metallic behavior of FeSe. 
Using an accurate single-particle description of the elec- 
tronic structure together with appropriate interaction 
parameters [l3|, and evaluating the influence of local 
Coulomb interactions within dynamical mean field the- 
ory (DMFT) [l^, we show that correlations are strong 
enough to give rise to the formation of Fe Zd local mo- 
ments, implying non-Fermi-liquid behavior, where elec- 
tronic states at the Fermi energy exhibit a finite life- 
time. Moreover, we demonstrate that these proper- 
ties are caused by a nearby doping-driven spin-freezing 
transition, i.e., Fermi- liquid behavior is restored towards 
larger electron doping, whereas hole doping reinforces 
bad-metallic properties. Using the same approach for 
FeAsLaO we show that Coulomb interactions are too 
weak to cause spin freezing, so that this system merely 
exhibits moderate effective mass enhancement [lij . 

Thus, FeSe appears to be a material which exhibits 
bad-metallicity induced by a spin freezing transition, a 
mechanism recently identify by Werner et al. 17| in a 
three-band model. For FeAsLaO, Haule et al. [l8| showed 
that large Coulomb interactions lead to the formation of 



local moments and large scattering rates, while Ishida 
and Liebsch p,9j discussed the spin freezing transition 
as a function of Coulomb interaction and doping. Bad- 
metallic behavior in FeSe was also found in DMFT stud- 
ies by Craco et al. [ii^l and Aichhorn et al. [2l|. The 
spin-freezing origin of this behavior, however, was not 
investigated in these works. 

To account for local Coulomb interactions among Fe 3(i 
electrons we use exact diagonalization [22J (ED) which 
we have recently extended to five orbitals [l^. Dis- 
cretization of the lattice surrounding the local impurity is 
achieved by using ten bath levels, which yields excellent 
projections of the lattice Green's function onto the clus- 
ter consisting of impurity plus bath. Because of the very 
large size of the Hilbert space of this 15 level system (the 
largest spin sector has dimension 40 x 10^) the spacing 
of excited states is very small so that finite size errors 
are greatly reduced. Moreover, ED has the advantage of 
allowing for rotationally invariant Hund exchange. This 
is of crucial importance for the spin freezing transition 
since omission of spin-flip and pair-exchange interactions 
leads to a significant shift of the Fermi-liquid to non- 
Fermi-liquid phase boundary. For computational reasons 
we restrict ourselves to the lowest excited states, which 
are relevant near the T — limit. Only the paramag- 
netic phase is considered. Further calculational details 
can be found in Ref. [l9| . 

The electronic properties of FeSe are formulated in 
terms of the effective low-energy model recently derived 
by Miyake et al. [IJ]. In this scheme, standard band 



structure calculations within the local density approxi- 
mation (LDA) were carried out, analogous to previous 
work in Refs. [2^, 24 1. From these results, maximally 



localized Wannier functions were constructed for the Fe 
3(i bands, following the procedure discussed in Ref. [25| . 
The inter-site transfer integrals are then derived from the 
matrix elements of the Kohn-Sham Hamiltonian -ff (k) in 
the basis of these Wannier functions. Finally, the con- 
strained random-phase-approximation (RPA) developed 
by Aryasetiawan et al. was used to determine the 
effective Coulomb and exchange interaction parameters, 
which account for screening via Se Ap orbitals. Because 
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FIG. 1: (Color online) (a) Low-energy scattering rates and (b) 
spin-spin correlations of FeSe as functions of 3d occupancy. 
Near n = 6, the non-Fermi-liquid (NFL) behavior 7^ > is 
caused by the formation of local moments. Increasing electron 
doping restores Fermi-liquid (FL) properties, whereas hole 
doping enhances bad-metallic behavior. 



of the planar structure of FeSe, these screening processes 
affect the various 3d orbitals differently so that the in- 
teraction parameters exhibit appreciable orbital depen- 
dence. According to these results, the key feature of FeSe 
is that Coulomb interactions are typically 50 % larger 
than in FeAsLaO. This follows from the larger z spac- 
ing and concomitant smaller spread of Wannier orbitals, 
and from the reduced number of screening channels in 
FeSe. The one-electron transfer integrals tmn(R) and in- 
teraction matrices Umn and j7m„ for FeSe are given in 
Tables VII and VIII of Ref. [l|], respectively. The av- 
erage intra-orbital Coulomb interaction is (U) — 4.2 eV 
and the average Hund exchange is (J) = 0.5 eV. The x, y 
axes point along Fe second-neighbor directions. 

Fig. 1(a) shows the low-energy limit of the 3d self- 
energy components, 7™ = — Iml]m(itj„ — > 0), as a func- 
tion of Fe 3d occupancy. The one-electron transfer inte- 
grals and interaction matrix elements are kept fixed at 
the values for n = 6. Ordinary Fermi liquid behavior is 
found for n > 6.5, but increasing low-energy scattering 
rates are obtained at lower electron doping, except for 
the dxy orbital. Even stronger scattering occurs on the 
hole doping side, n < 6. The reason for the different be- 
havior of jxy is that pxy{uj) is about 0.5 eV wider than 
the other density of states components 

To identify the origin of this change from Fermi- 
liquid to non-Fermi-liquid behavior we have evaluated the 
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FIG. 2: (Color online) Fe 3d orbital occupancies as functions 
of chemical potential. Uav denotes the average occupancy. 
Nominal occupancy n = 6 is reached near /i = 17 eV, while 
at /X = 15 eV all bands are approximately half-filled. 



spin-spin correlation function C„i(t) — {Smz{T)SmziO)), 
where r denotes imaginary time. In the case of a Fermi- 
liquid, these functions decay with t, so that Cm(r — (3/2) 
is very small (/3 = 1/T = 100 eV^^). The susceptibili- 
ties Xm ~ dr {SmziT)SmziQ)) then are Pauli-like, i.e., 
independent of temperature. Fig. 1(b) shows that the 
results for n > 6.5 are consistent with this behavior. At 
smaller electron doping, the Cm=^xy{T) components reach 
increasingly larger constant values near r = /3/2. The 
corresponding susceptibilities Xm then become propor- 
tional to 1/T, as expected for Curie- Weiss behavior. 

A similar spin-freezing transition was recently found 
by Werner et al. for a degenerate three-band model 
near n = 2, i.e., at one electron away from half-filling. 
Using continuous-time quantum Monte Carlo DMFT, the 
paramagnetic phase at moderate U was shown to exhibit 
Fermi-liquid properties at small n. For n > 1.5, an in- 
coherent metallic phase appears, where the self-energy 
exhibits a finite onset at lo ~ due to the formation of 
local moments. Because of particle-hole symmetry, the 
same Fermi-liquid to non- Fermi-liquid transition appears 
on the electron-doped side, i.e., when reducing the occu- 
pancy below n w 4.5. 

Fig. 2 shows the 3c? orbital occupancies as a function 
of doping. At n = 6, these values differ significantly 
from their uncorrelated values: nxz^yz,x^-y^.xy,z^ — 
(0.54, 0.54, 0.56, 0.60, 0.76), indicating strong interorbital 
charge transfer induced by Coulomb interactions. In par- 
ticular, rixy reaches nearly unity, while the other Tijji are 
close to half-filling. Despite this orbital polarization, the 
spectral distributions (not shown here) reveal consider- 
able intensity near Ep, suggesting that the system is not 
in a Mott phase, where one of the subbands is filled and 
the others are split into lower and upper Hubbard bands. 
These kinds of partial Mott transitions, with certain sub- 
bands empty or full, are found in several three-band t2g 
materials, suchas LaTiOs [13, [ll], V2O3 [H IH, and 
Ca2Ru04 
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To investigate the relationship between orbital polar- 
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FIG. 3: (Color online) Low-energy scattering rates and spin- 
spin correlations derived within ED/DMFT for a degenerate 
five-band model with Bethe lattice density of states with W — 
4 eV. For [7 = 4 eV, J — 0.75 eV, the spin-freezing transition 
occurs at n = 6.2, with Fermi-liquid behavior at n > 6.2 
and increasing bad-metallicity for n < 6.2. For (7 = 4 eV, 
J — 0.5 eV (not shown), the transition occurs at n = 6.1. 



ization and local moment formation, we have calculated 
the effect of Coulomb correlations in a fully degenerate 
five-band model consisting of identical densities of states 
for a Bethe lattice, with band width and Coulomb in- 
teractions of similar magnitude as in FeSe. As shown 
in Fig. 3, this model system exhibits the same kind of 
spin-freezing transition as discussed above for FeSe. For 
n = 6, the low-energy scattering rate indicates non- 
Fermi-liquid behavior, associated with the formation of 
local moments. Fermi-liquid behavior is recovered at 
electron doping n > 6.2, whereas hole doping strengthens 
bad-metallic properties. The important conclusion from 
this result is that orbital polarization is not a prerequi- 
site for the spin-freezing transition. Thus, this picture 
differs from the orbital-s elective , itinerant-localized sce- 
narios proposed in Refs. 



To understand better to what extent the results for 
FeSe shown in Fig. 1 depend on the details of the 
single-particle Hamiltonian, we have carried out anal- 
ogous ED/DMFT calculations by replacing i/(k) with 
the tight-bin ding Hamiltonian derived by Graser et al. 
for FeAsLaO [33| . Since the 3(i density of states compo- 
nents of this compound are qualitatively similar to those 
of FeSe, there ought to be also a spin-freezing transition. 
Indeed, if we retain the FeSe interaction parameters with 
an average Coulomb energy of 4.2 eV, the overall be- 
havior of the low-energy scattering rates 7m and spin- 
spin correlations C,„ are qualitatively similar to the ones 
shown in Fig. 1, except for a less pronounced orbital po- 
larization. For instance, at n = 6, n 



xy 



0.77 rather 

than 0.95. Nevertheless, the Fermi-liquid to non-Fcrmi- 
liquid transition also occurs near n — 6.5, suggesting that 
this phenomenon is remarkably robust, i.e., insensitive to 
the details of the one-electron properties. These results 
further support the key point of the degenerate-band ex- 
ample presented in Fig. 3, namely, that the spin freezing 
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FIG. 4: (Color online) (a) Low-energy scattering rates and 
(b) spin-spin correlations of FeAsLaO as functions of 3ci oc- 
cupancy. The Fermi-liquid to non-Fermi-liquid transition 
caused by local moment formation occurs close to n = 6. Elec- 
tron doping maintains Fermi-liquid properties, whiles hole 
doping leads to bad-metallic behavior. 



transition is not driven by orbital polarization. 
According to the work by Miyake et al 



1J|, the cru- 
cial difference between FeSe and FeAsLaO is that, be- 
cause of more efficient dp screening and more extended 
Wannier orbitals, local Coulomb interactions in the lat- 
ter system are considerably smaller, with {U) « 2.7 eV 
and (J) 0.4 eV. In previous ED/DMFT calculations 
for FeAsLaO [l^] based on orbital independent Coulomb 
and exchange energies we showed that spin freezing may 
occur at about C/ = 3 eV for Hund exchange J — 0.75 eV. 
Using the recently published orbital dependent Umn and 
J„in matrices derived within constrained RPA 



14 



we 

are now able to make a more accurate prediction of the 
spin freezing transition in FeAsLaO. The transfer inte- 
grals and interaction parameters are again held fixed at 
their values for n — Q. As shown in Fig. 4, the tran- 
sition now is located almost exactly at n = 6. Thus, 
in contrast to FeSe, Fermi-liquid properties in FeAsLaO 
prevail and the 3d bands exhibit only moderate effec- 
tive mass enhancement, m* « 2... 3 [16j. Nonetheless, 
hole doping n < 6 gives rise to local moment formation 
and bad- metallic behavior, with a finite lifetime at Ep, 
while electron doping stabilizes the Fermi-liquid proper- 
ties. (Note the weaker orbital polarization compared to 
the results for FeSe in Fig. 1. Also, because of orbital de- 
pendent Umn, Jmm the magnitudes of 7m and C'miP/^) 
differ from those for orbital independent U, J in Ref. [ISj.) 
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FIG. 5: (Color online) Schematic phase diagram for FeSe and 
LaFeAsO. Solid curve: spin freezing transition indicating the 
boundary between Fermi-liquid and non-Fermi-liquid phases. 
Vertical bar at half-filling: Mott insulating phase. 



On the basis of these results we arrive at the phase 
diagram depicted in Fig. 5. FeSe is located well inside 
the non-Fermi-liquid phase, whereas FeAsLaO lies on the 
Fermi-liquid side of the spin-freezing transition. Both 
systems exhibit 'parent' Mott phases in the limit n = 5. 

We finally discuss the role of spin-flip and pair- 
exchange interactions. In Ref. 19] it was shown that 
the omission of these terms in the many-body impurity 
treatment leads to a significant shift of the Fermi-liquid 
to non-Fermi-liquid phase boundary to ^ 1 eV smaller 
(TJ). Thus, for FeSe bad-metallic behavior would become 
even stronger. Since FeAsLaO, however, lies very close 
to the phase boundary, this approximation has severe 
consequences. Instead of a Fermi-liquid characterized by 
moderate effective mass enhancement, this compound is 
then also bad-metallic, in conflict with experiments. 

The present work is based on a consistent combination 
of single-particle Hamiltonian, constrained RPA interac- 
tion parameters, and DMFT many-electron calculations 
within a 5 X 5 d electron basis. It would be interesting to 
compare the spin freezing transition obtained within this 
scheme to analogous formulations within a pd electron 
basis, 



21|. 



such as the one used in Ref. 
In summary, we have evaluated the effect of local cor- 
relations in FeSe, using accurate single-particle proper- 
ties and constrained RPA results for the orbital depen- 
dent Coulomb and exchange matrices, combined with 
ED/DMFT. The finite scattering rates derived from the 
3d components of the self-energy are shown to be linked 
to the formation of local moments. Thus, FeSe is bad- 
metallic, in agreement with experimental findings. Spin 
freezing ceases for large electron doping, giving rise to 
Fermi-liquid behavior. In contrast, hole doping enhances 
bad-metallicity. Applying the same scheme to FeAsLaO 
we find that, as a result of the more efficiently screened 
Coulomb interactions, local moment formation is con- 
fined to hole doping n < 6. Thus, FeAsLaO and FeSe 



seem to be located on opposite sides of the Fermi-liquid 
to non-Fermi-liquid spin freezing transition. 
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